DGP black holes on the brane 
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We find an exact solution on the brane for static black hole in the DGP model. In the appropriate 
limit we recover the 2 known solutions, the Schwarzschild and the Reissner-Nordstrom solutions 
with tidal charge. The solution has 2 branches, which correspond asymptotically to de-Sitter or flat 
Universe. Also we study the linear stability of the solutions. We find that the Regge- Wheeler and 
the Zerilli potential are positive. Finally the dipole perturbation is derived which corresponds to a 
linearization of a rotating black hole solution on the brane. 
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Extra dimensions provide an approach to modify grav- 
ity without abandoning the form of the action proposed 
in the Einstein's general relativity. From a phenomeno- 
logical point of view we can avoid constraints com- 
ing from standard model observations, by considering a 
brane world scenario, that is, we are living in a hypersur- 
face (3-D) in a higher dimensional spacetime. From the 
theoretical point of view, string theory predicts a bound- 
ary layer, a brane, on which edges of open strings stand 
[H . The possibility that we may be living in a brane gen- 
erates many questions as how gravity looks like. Also in 
an attempt to solve the much debated hierarchy problem, 
various problems are studied, but also in order to under- 
stand the cosmology, such as inflation and Dark Energy. 
In this contribution to consequences of the brane-world 
in 4-D, we study one of the most famous model, DGP 
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We know that the physics of black hole and gravita- 
tional collapse is complicated, especially because of the 
matter localized to the brane, while the gravitational field 
can access the extra dimension, and also because of the 
nonlocal effects of the bulk into the surface. Part of the 
problem would be to find a solution on the brane. We 
can not necessarily embed this solution into a bulk but 
some information of the global solution can be under- 
stood, some intuition can be developed. The solution 
can be smoothly continued into the bulk via the ADM 
formalism, where the solution on the brane can be con- 
sidered as an initial data. At least a local solution of the 
bulk exist, even if the global solution is not guaranteed. 

The model is defined as an empty five-dimensional 
space (not necessarily Minkowsky) and all the energy- 
momentum is localised on the four-dimensional brane. 
The theory is described by the following action in the 
vacuum 



S 




(1) 



where (g, h) are respectively the metric of the bulk and 
the brane and R the intrinsic curvature in 5-D and 4-D 
respectively. The Gibbons-Hawking term is implicit in 

©• 

Variation of the action gives the following equations of 
motion 



Bulk equation: iJ„„ = 0, 



Brane equation: Guv 
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(2) 
(3) 



where r c = MhJ2M%^ is the crossover scale that gov- 
erns the transition between four-dimensional behaviour 
and five-dimensional behaviour. 

Following i and i we can rewrite an equation on 
the brane with the metric h only. For this we define the 
tensor H 



K 1 



which gives on the brane 



(4) 



(5) 



and L is solution of the following algebraic equation 
(Gauss equation) 



T 2 ^ 

T a T - —h 4- h 



-E^, (6) 



where L is the trace of the tensor and E is the electric 
part of the Weyl tensor. 

In the following, we will focus on static spherically so- 
lutions in the vacuum of the form 



ds 2 



-A(r)di 2 



dr 2 



(7) 



Therefore we can decompose irreducibly the tensor E 
with respect to a 4- velocity field m m [H, [f| 



E^u = p{r) (u^u u + + P(r) [r^r v - jjjg^) , (8) 

where q^ u = h^ u + u^u v projects orthogonal to the 
timelike vector. Here (p, P) are respectively an effective 
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energy density and anisotropic stress on the brane arising 
from the 5-D gravitational field. 

It is easy to sec from (|5|7j) that L is diagonal, hence 

we write U l u = di&g(^L , L l7 L 2 , L^. 

From the equations (|6I8[) we have the following alge- 
braic equations 



we can close the system of equations and get an equa- 
tion for v 



d^-3)(^-3)=0, 



(16) 



L 1 = ±dL 2 --(2p + P), 



Lo — — ± 



r 2 



(9) 
(10) 



where the equation L 2 = L-$ comes from the Eq.([5]). 
Also from the brane equation ([5]), we have 



(ii) 



Hereafter, we will assume that radial photons should 
experience no acceleration, the velocity of light in the 
radial direction should remain constant. Therefore we 
have Q A = B which implies L\ = L$, hence from © we 
have 2p+ P = 0. This constraint between the density and 
the pressure is the same than in the absence of induced 
curvature term Q. 

Finally, we have from ^ 



4L 2 + 2P±8L Jl 2 + P 



(12) 



Following Q, we define v = 2 ± V ' L\ + P/L . Then, 
it is straightforward to check that 



T 2 - 



P = 



2r 2 c (v 2 -3)' 
3 (v-l)(v 



3) 



2r 2 



(13) 
(14) 



which gives 



|^_73|(^3+D/2 

\v + VSl^- 1 )/ 2 ^ -3| : 



(17) 



where Q is an integration constant. 

The coefficient r 2 is fixed in order to recover the result 
P ~ 1/r 4 Q in the limit r c — > 0, which would corre- 
sponds to the Reissner-Nordstrom solution on the brane 
with tidal charge. 

The eq.tfTT]) gives v(r) which from ([13)) gives Lo(r), 
therefore we can solve the final equation ([5]) 



"B'(r) +B(r) - 1 



2r 2 . 



± 



(18) 



where the sign ± is because of the sign of Lq . 

We have found 3 different solutions depending on the 
range of v. In the first solution we have v < — \/3, the 
second solution corresponds to y/3 < v < 3 and the last 
to v > 3. Accordingly the range for r will be respectively 
r > \JQr c , r > and r > y/Qr c . The second and the 
third solutions are identical except the range for r, hence 
we keep only the second solution which covers the full 
spacetime (brane). The first solution do not cover the 
full spacetime r > \JQr c and therefore can't describe 
a black hole. This solution will not be studied in this 
paper. Hence we have 2 branches of the solution 



We see that we need v 2 > 3. 

The only undetermined function is v, all the other 
quantities as (p, P, K^ v ) are related to v. 
Considering now the Bianchi identity 



V M Lt + U/ V L = 0, 



(15) 



A = B = 1 



r 2r 2 . 2r 2 J^ ' 



(19) 



where m is an integration constant, v is solution of the 
algebraic equation p7| and / can be written in terms of 
Gauss's hypergeometric function, 



f(v) = v/6 



Vv 1 



4V2"(n/3-1)(v-3) 2 
5(«- V3)V 2 (v + V3) x /2 

r 



iFi 



3(3 



(20) 



/ is a concave down monotonically increasing function of r. It is negative for r < 0.78\/Qr c and positive other- 
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wise and lim r _ ) . 0O / = 1. 

The horizon structure of the black hole on the brane 
depends on the branch considered. The negative branch 
(negative sign in the equation ([TT)]0 has a black hole hori- 
zon and a cosmological horizon because of the de-Sitter 
structure at large distances, while the positive branch, 
which is asymptotically flat, has a single horizon. 

It is interesting to study the asymptotic behaviour of 
this solution. We have at large distances 



f(r) = 1 - 



2q 2 



2q* 



6^6 



5r 8 



9r 12 



+ OM, (21) 



where we have redefined the integration constant q = 

V / 372(3 - V3) (v ^~ 1)/4 (3 + V3)~ (V ^ +1)/4 Q- 

We see that only the positive branch has a smooth 
limit as r c — > (Randall-Sundrum model limit) and as 
such we will refer to as the RS branch. In contrast the 
negative branch is not smooth as r c — > and represents a 
distinct new feature of DGP, the DGP branch also known 
as the self-accelerating branch. The RS branch converges 
to A = 1 — 2m /r — q 2 /r 2 , it is not a Reissner-Nordstrom 
spacetime. In fact the tidal charge has always the same 
sign and it is physically more natural for a brane solution 
Q, the tidal charge strengthens the gravitational field. 
This is why our solution do not have a Cauchy horizon, 
even in the limit r c — > 0. Also as r c — > oo, we recover the 
Schwarzschild solution (1— 2m/r) for the 2 branches, the 
solution is the same than in Einstein theory, there is no 
vDVZ IH discontinuity. The continuity of the theory 
is restored because the nonlinear effects were taken into 
account while we would conclude to a discontinuity if we 
use the linearized solution at large distances (f2"Tj) . 
At small distances 



f(r) 



a 







W /3(l + v^)/8 ^(3V3-5)/8' 

where (a, /3, 7) arc 3 positive constants 



(22) 



3 (2^ - 3) 



\(V3-3) 



r (I) r (I (373-1)) 



a/45 + 26^ T (| (3 



V3)) 



0.76, 



lll\/3 



1.05, 



3^3 (2V3 - 3) " T (f) T (i (3^3 - 1)) 



5\/2r 



V3)) 



0.77, 



and W is the Lambert function at the point 

3 V3/2 4 l-V3( 2 -^)( 



IVc 



The dominant contribution is given by the first term 



f(r) 



\r c qfl 2 



(23) 



it* 




FIG. 1: Existence of the black hole for the DGP branch: The 
blue part represents the range of the parameters (m, q) for 
which the black hole has 2 horizons. In the grey region, the 
metric is always negative and the white region corresponds to 
a black hole with only a cosmological horizon, because of the 
negativity of the effective mass, hence inside the cosmological 
horizon the solution would be a naked singularity. 



with 



2^2- 



r(f)r(|( 3V3-i)) 

3 + V3)) 



3 3 / 8 r 



3.11. 



We see that even in the case of a massless black hole, 
we have a "mass-term" because of the fifth dimension. 
Hence we recover a standard result; even for a massless 
black hole, the behaviour of the solution is 1/r at small 
distances and 1/r 2 at large distances. 

In order to keep the effective mass positive, in the DGP 
branch, we impose q < |m, where (m = m/ ^lqr c , q = 
q/r c ). The existence of the black hole is constrained in 
the Fig.([T]). From which we see that for a fixed parameter 
q, the mass of the black hole has an upper bound but also 
a lower bound. We would have a naked singularity for 
lightest black holes, this can be seen as an instability of 
the branch. 

The RS branch is much simpler, we have a black hole 
for all positive parameters (fh,q), also these parameters 
play the same role, they increase the position of the hori- 
zon and hence its entropy. At large distances, the New- 
tonian potential is dominant. Depending on the param- 
eters, the situation can be the same for all distances, the 
solution will be very close to the Schwarzschild solution, 
except for large values of q where the mass of the black 
hole will be renormalized at small distances. But in the 
case of the DGP branch, we do not have the same be- 
haviour at large and small distances, hence we have a 
new distance scale r+ dubbed the Vainshtcin radius 
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r* - K) 1 / 3 . (24) 

As we said previously, the mass of the black hole is 
bounded from bellow. The existence of the horizon is 
constrained by m > g 2 r^/r* ~ g 2 10 28 A'/ Q if we assume 
the Vainshtein radius of the order galaxy scale and the 
crossover scale of the order Hubble scale. A stellar black 
hole exists if q < 10 -14 . Otherwise it will be a naked 
singularity. 

The extrinsic curvature and therefore the curvature 
constant can be easily derived from (|13I14|) 

i?=A( 2T V6^=L), (25) 
ri V Vf 2 -3/ 

which is singular at r = 0. Also we can see that 
R > for the DGB-branch and converges to 12/r 2 , 
while we have R < for the RS-branch and it goes 
to zero at infinity. A non-vanishing curvature outside 
the source leads usually to screening mechanism and 
we have shown previously the absence of the vDVZ 
discontinuity. On the physical stability of the solutions, 
it is interesting to study the violation of the energy 
conditions if we consider the tensor K^ v — Kh^ v as a 
source term, and the positivity of the gravitational mass 
of this spacetime. These particular problems should be 
addressed separately. 

In order to study the linear stability of this solution, 
we follow the Regge- Wheeler formalism [13, [ll[ and we 
decompose the metric perturbations according to their 
transformation properties under two-dimensional rota- 
tions. They are classified depending on the transforma- 
tion properties under parity, namely odd (axial) and even 
(polar). Using the Regge- Wheeler, and Zerilli gauge, one 
obtains two distinct perturbations : odd and even per- 
turbations. 

For i > 1 , the equation of perturbations takes the form 

^~^ + V^ = S, (26) 

where S is the perturbation of the source term, r* is 
the tortoise coordinate defined as dr* 2 = dr 2 /A, ip is a 
function of the metric perturbations and V is the Regge- 
Wheeler potential or the Zerilli potential in the respective 
case. 

We have for the Regge- Wheeler potential 

_ A (X + 2A- rA') Ylpr 2 A 

RW r 2 6r 2 + r 2 \{r 2 A)" - 2] ' 1 ' 

W- 1 + 3,4 3A / I 6 \ 

= A ^^ + 27l( lT V^)' (28) 

where A = (£ - l)(i + 2). 



The last form is useful in order to see that the potential 
is always positive for the 2 branches. The first term of 
the potential is the standard term in 4 dimensional and 
the last term comes from the 5 dimensional effects. We 
do not write here the Zerilli potential because it is much 
more complicated but as in General Relativity, the graph 
of the Zerilli potential is similar than the Regge- Wheeler 
potential for the same parameters. 

The positivity of the potential indicates the stability 
of the spacetime under linear perturbations fl2j for the 
2 branches. It is important to notice that we didn't con- 
sidered the source term in order to study the stability of 
the theory. In this case, the source term is much more 
complicated than in General Relativity. In fact it is not 
localized, the source is function of the electric part of the 
Weyl tensor which is a non local term. 

For the dipole perturbation t = 1, we can write 
ht(j> = P{r) sin 2 (#) because the time dependant term can 
be removed via the gauge freedom. The equation for j3 
takes the form 

F=1±S (29) 

where we have neglected the source term. The solution 
of this equation represents the spacetime around a slowly 
rotating black hole on the brane for the DGP model. In 
the case where r c — > oo which corresponds to the GR 
limit, we recover the standard result, the Kerr metric at 
the first order of perturbations. 

h tq > = -^sm 2 (6), (30) 

where J is identified to the angular momentum and we 
have gauged away the unphysical term proportional to 
r 2 . 

The Eq. (j2T)|) has a Hill type form, hence the evolution 
will depend on the sign of the RHS of the equation. In or- 
der to integrate the system, we assume a Kerr form of the 
metric at small distances. For the DGP branch, we have 
an angular speed f2 = —/3(r)/r 2 which decreases with 
the radial distance but contrary to General Relativity, 
il = at a finite distance, smaller than the cosmologi- 
cal horizon. After this critical point, the frame-dragging 
is occurring in the opposite direction. In the case of the 
RS branch, the angular speed depends on the parameters 
(rh,q). For small q we recover Kerr solution. Once we 
increase this parameter, f2 goes to a finite value at infin- 
ity, and finally when q > rh the angular speed increases 
after decreasing and converge also to a finite value. 

In conclusion, we have derived an exact black hole so- 
lution on the brane for the DGP model. This solution 
recovers the standard results at small and large distances 
but also covers the intermediate regime which was not 
known. The 2 branches depend on 3 parameters, the 
mass of the black hole, the tidal charge and the cross 
scale parameter. We have shown than if we do not con- 
sider the perturbations of the bulk (the source term), the 
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solutions are stable under linear perturbations. It would 
be also interesting to see how the horizon will be affected 
by the full solution and how the quasinormal modes of 
the black hole are modified. Finally we have shown from 
the dipole perturbation, an approximation of a possible 
rotating black hole on the brane. For the DGP branch, 
the metric g tt j, decreases to zero in a finite range of the 
radial distance, while for larger distances, the angular 
speed is opposite. Whereas g t <f> converges to a finite value 
at infinity for the RS branch. Contrary to General Rela- 
tivity, the dragging do not disappears completely at large 



distances. 
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